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SUMMARY 


The performance of Model Reference Adaptive Control (MRAC) is 
studied in nSerLa! simulations with the S3, 

the effects of differences between the pi ant and the refe 

MRAC is applied to two structural sys ems £ 3 Performa nce indices 

KoflI PDn f hp reference model and the actual piano. 

relating^”"? oontrojef fort and response oharaoteristios are monitored 
S Sde? to determine what effects small errors have on the control 
effort and performance of the two systems. It is shown that 
reasonable amounts of error in the reference model can cause dramatic 
increases in both the control effort and response magnitude (as 
measured by energy integrals) of the plant 


INTRODUCTION 


nurina the past decade, researchers have shown much interest in 
control and wSntlllcation of large flexible structures with emphasis 
on Large Space Structure (LSS). Furthermore, our inability Jr®.”®' 1 ® 
these large structural systems accurately has generated extensive 
research into adaptive controllers capable of “ a jf^ n ^ n ^ lllty in 
the face of large structural uncertainties as well as changi g 
structural characteristics . However, most of this research has been 

strictly theoretical in nature (eg . refs 1-10, 

far ‘behind?" in'ldai^ of mLt' ‘theoret ical research has 

been on designing stable adaptive controllers with little or no 
concern for the issue of control effort. 

While it is possible to design an adaptive controller that will 
stabilize a structure even if we have a very poor model, the control 
effort may be very high. The objective of the present paper ^ to f 
study the Y correlation between the control effort and the fidelity o 
Se strSrtSS model . Specifically, the first step ^ “ le,o„str a te 
that the effort associated with an adaptive control system 
sensitive to knowledge of the structure. For this purpose the; P°P 
Model Reference Adaptive Control (MRAC) method 3 eleoted and two 

examples were studied in detail. In this paper, we monitor four 
performance indices: Maximum Control Force, Quadratic on ro , 
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Kinetiic: Energy and Potential Energy. These performance indices 
allow us to evaluate the effects of errors in the theoretical model 
Numerical simulations were used to see how each performance index 
changed when errors were introduced into the system. Section II 
summarizes the MRAC algorithm, section III shows how a simply 
supported beam can be sensitive to the choice of Reference Models 
sections presents the sensitivity of a more co»pUc”ed sSuc&e 
and section V provides concluding remarks. 


MODEL REFERENCE ADAPTIVE CONTROL 


Adaptive controllers generally fall into two classifications 
direct and indirect. The basic difference between the two 
lassifications is system identification. Indirect adaptive methods 
• 9 1°' 1 3 ~ 14 ) require system identification before the 
% 9a i nS thS controller car > be updated, whereas direct 
^ %i refs * 1-8 ' 11-12 ) do not use system identification. MRAC is 

adapt ivelv lunt P°P ular di f? Ct methods (refs. 1-7). MRAC methods 
_ P tlvely the contr oller gains, forcing the actual system to 

S ° me ld ^ al £® ference model. Because this reference model can 
be .° f . 1 ? wer order than a typical model of the actual system this 
mSd^o i! V ; ry f^ctive for applications to LSS, wheJe struct Jral 
anfinnJJon 6 ° f hlgh order and require truncation for use with 

system (ref.^) igure 1 shows a block diagram of a generalized MRAC 


PROBLEM FORMULATION 


spaced Srm? S ' ° T controlled P lant can be represented in standard state 


X p (t) - ApXp(t) + BpU p (t) 


(la) 


Y p(t) - C p Xp(t) (lb) 

where X p e R Np , U p e r m , y e r m and a , B , C are of appropriate 
dimensions. it is assumed that (A , B ) P is P controllable (A c ] is 

and that thS " U " ber ot in P uts (») equal to tfie n’umber 

s t abl e reference model which specifies the desired performance 
of the plant is also described by a state space representation, 
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X m (t) - VUt) + B m U m (t) 


(2a) 


c m x m (t) 


(2b) 


where X m € R N \ U m 
dimensions. For 
must be true 


« rM ' Y m 6 
practical 


R h and A , B m , C are of appropriate _ 
application to LSS the following condition 


N„ > > N„ 


(3) 


To aid in measuring how close the actual plant is to the reference 
model, the output error between the plant and the reference model 

defined as 


e y (t) - Y m (t) 


Y D (t) 


(4) 


Since the output error tells us how close the actual plant is to the 
desired performance of the reference model, the objective of any 
adaptive update scheme is to design a control input which forces the 
output error to zero asymptotically. 


ADAPTIVE LAW 

The adaptive mechanism used in this paper is based on the work of 
Sobel , Kaufman and Mabius (ref. 1) and has no provisions for the 
destabilizing effects of noise. However, it has found important 
application to LSS. The control input is written as, 

U p ( t ) - K ( t) r (t) ( 5 ) 


where 



(6a) 


K ( t ) - [K e (t) ,K x (t) ,Ky(t)] 

The adaptive gain K(t) is calculated as the sum of a proportional 
component K pr (t) and an integral component K,(t) so that 

K (t ) - K pr ( t ) + K[ ( t) 

The adaptive laws for K pr (t) and K t (t) are given as 

K pr (t) - e,r'T- 
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K,(t) - e y r T T 


(9) 


w here T and T are time invariant weighting matrices of appropriate 
dimension chosen by the designer. Sufficient conditions for global 
stability are presented in (ref 1-3, 11-12) and will only be 
summarized here. 

1. T & T* > 0 

2. there exist P = P T > 0 and Q = Q T > 0 such that 


PA p + ApP - -Q 

Condition 1 is met simply by choosing appropriate matrices (i.e., the 
identity matrix) . Condition 2 is equivalent to the assumption that 
the open-loop plant transfer function matrix 

Z(s) = C p (s I - A p ) ’B p (10) 


is strictly positive real. This condition is met for any LSS having 
small but non-zero inherent damping and colocated sensors and 
actuators . 


CONTROL EFFORT 

In order to assess the added implementation costs of MRAC in 
systems where reasonable amounts of error would occur, we have adopted 
the following procedure. The first step is to choose a linear system 
to represent the accual physical system. Next, we create a reference 
model which specifies the desired performance and has some measurable 
amount of error. Previous examples, see (ref. 2), have chosen the 
reference model to be a reduced model of the actual plant with the 
same frequencies and mode shapes plus extra damping. While this would 
be the ideal situation, it is not probable that we would have an exact 
theoretical model. For this reason we have intentionally introduced 
errors between our reference models and the actual plant model. To 
aid in quantifying the increased effort due to the errors, we 
calculate the following performance indices: 

the maximum control force required by each actuator. 
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uf x = Max ( | U,(t) I ) 0 < t < t final , i - 1, M 


( 11 ) 


the quadratic control effort, 

Utotal - JU T Udt 


(12) 


the integral of the potential energy of the system. 


PE 


X T KXdt 


(13) 


and the integral of the kinetic energy 

KE - J X T MXdt ( 14 ) 

where the first two performance indices measure the control effort and 
the second two provide information about system response 
characteristics. These performance indices allow us to see how 
increments of error affect the cost and performance of the system. 


SIMPLY SUPPORTED BEAM EXAMPLE 

The first example is a simply supported beam with a variable 
concentrated mass at the mid-span and a velocity sensor and force 
actuator colocated at one-sixth span (see figure 2) . is simp ® 
structure is similar to a structure used by Bar-Kana, Kaufman & Balas 
(ref. 2) for demonstrating the MRAC method. The only difference 
between the present structure and the structure of reference 2 is e 
variable concentrated mass. The variable concentrated mass at the 
mid-span was used to create error in the system due to unknown mass 
characteristics. The concentrated mass was varied between 0 2 0 -s of 
the mass of the beam, with zero mass corresponding to an exact 
reference model. It should be noted that the reference model was 
held constant while the plant model was varied to match changes in t 
concentrated mass. 

The beam was modeled with 12 beam finite elements with a 
displacement and rotational Degree of Freedom (DOF) at each nodal 
point. The coupled equations of motion are written in standard form 

as 


Mq + Cq + Kq - F 


( 15 ] 


where M and K are the mass and stiffness matrices respectively, and C 
is the damping matrix calculated from assumed inherent damping ratios 
£.. Using modal analysis the equations are transformed from a se o 
coupled equations in physical coordinates to a set of 
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coupled equations in physical coordinates to a set of uncoupled 
equations in modal coordinates 

X + 2CwX + AX - B°F (16) 


where 


A - diag [c^, <£, . . . c^ 2 ] 


the ©j's are the undamped natural frequencies 

C “ diag [ C i W| t C z^z / • • • / (12 ^12 J 


B “ [^5,1 '^5,2' • • • »^5,12] T 


and the 0c i ' s are the fifth element of each eigenvector (the sensor 
and actuator are at the 5 th DOF) normalized so 

# t M$ - I 


Equation 16 is rewritten in state space form as 


*p 


-2Cp^p -Ap 

*P 

+ 

Bp 

*p. 


I 0 

X P. 


— 

O 


(17a) 



(17b) 


where the subscript p denotes the equations apply to the plant. The 
reference model takes the same form, 


X m 


2 C m t 4wn “An 

K 


b° 





+ 

m 



1 0 . 

x* 


0 . 


(18a) 



(18b) 


where the subscript m applies to the reference model. 

For the purpose of numerical simulations we must reduce the size 
of the actual plant. In this example (as in ref. 2) we consider only 
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the first three modes of the actual plant and choose a reference model 
that includes only 2 modes. Damping ratios in the plant are assumed 
to be 0.01 while the desired damping ratios of the reference model are 
set at 0.05. All other parameters (length, El, etc.) are set to 1.0 
for convenience in calculations. In the present study we consider 
only initial condition responses. The first three modal states were 
initially set to 1.0, while all others were set to zero. From figures 
3-4 it can be seen that the controller does an excellent job of 
forcing the actual plant to follow the reference model. However, from 
Table 1 it can be seen that the addition of the concentrated mass, 
i.e. errors between the actual plant and the reference model, can 
produce very large increases in the maximum control force and control 
effort needed for the controller to function. For example, a 
concentrated mass weighing 20% of the beam weight causes a factor of 
six increase in the quadratic control effort. This large increase in 
control effort demonstrates a need to find a method for choosing a 
good reference model . 


SLEWING GRID EXAMPLE 


The Virginia Tech slewing grid laboratory structure shown in 
figure 5 is a more complex example. The slewing grid was designed to 
have characteristics of LSS, namely closely spaced modes, low 
vibration frequencies, and low inherent damping. Three pairs of 
velocity sensors and force actuators are colocated at joints 3,4 &5. 
The slewing grid was designed to include a zero frequency rigid body 
rotation mode about the shaft, but this has never been realized 
because of bearing friction and slight misalignments of the rotational 
shaft. Although the geometry of the structure is symmetric about a 
horizontal line through joint 3, the vibration mode shapes are not 
similarly symmetric because the structure's weight causes asymmetric 
member gravity loading and therefore asymmetric stiffness 
distribution. It was considered desirable in the design phase to have 
at least one beam member in substantial compression relative to its 
buckling load, both to reduce the overall structural stiffness and to 
permit the possibility of nonlinear response. The lower horizontal 
member carries the largest compressive load, being compressed to about 
70% of its Euler (pin-ended) buckling load. Great effort has been 
taken to accurately predict the loads in each member of the structure. 
However, each joint is held in place with a nut and bolt assembly and 
the process of tightening these bolts induces forces which we have 
been unable to determine accurately. Therefore our current Finite 
Element Model (FEM) only takes gravity forces into account. The 
rotational shaft was modeled by 8 beam finite elements with a 
displacement and rotation DOF at each node. Each of the 5 members of 
the structure is modeled with 4 finite elements which include a 
transverse displacement, an in-plane rotation, and an out of plane 
rotation at each node. The complete FEM has 72 DOF and the coupled 
equations of motions can be written 
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(19) 


MX + CX + (K + G) X = F 

where G is the geometric stiffness matrix. To make the problem more 
manageable we created a reduced eleventh order model using the Guyan 
Reduction (ref. 15) . The linear eguations for the slewing grid can 
be written (in physical coordinates) 


where 


X„ 


x n 


-M C 

u p '-p 


- M pX 


0 

k 


* 

b! o 

p 

L p 



. pj 

Kp ♦ ■ 

G P 


[EL 


u„ 


(20a) 


(20b) 


( 2 1 ; 


M p , C p , K p , and G p are the reduced mass, damping, stiffness and 
geometric stiffness matrices, and B p is (22 X 3) matrix with only 3 
non-zero elements for mapping the control inputs to the proper DOF at 
joints 3,4,5. 

The accuracy of the frequencies and modes predicted by the FEM is 
not good. During the past two years, great pains have been taken to 
find a FEM which would accurately model the structure. However all 
the non-linearities in the structure, such as friction in the bearing, 
large gravity loading in the lower horizontal member, and the loads 
induced by tightening the bolts at each joint of the structure, have 
resulted in a modeling nightmare. Table 2 compares frequencies 
predicted by our best FEM to the experimental vibration frequencies, 
and figures 6-9 compare several experimental and theoretical mode 
shapes. The difficulty in accurately modeling this structure was the 
driving force behind the decision to apply adaptive control to the 
slewing grid. The challenges of modeling the slewing grid may be 
similar to those we will face when we begin to model LSS. 

In order to study the performance of MRAC for this case we first 
had to choose a model for simulating the actual plant. Our efforts to 
model the slewing grid as accurately as possible resulted in several 
FEM with varying degrees of accuracy. The most accurate model used 
experimental frequencies and mode shapes in a correction method 
proposed by Baruch (ref. 16) to force the theoretical model to have 
exact experimental frequencies. The least accurate model was the 
standard FEM with no corrections. With this in mind we chose the most 
accurate FEM to simulate the plant and a linear combination of the 


426 



The 


most accurate and least accurate model as the reference model, 
reference model is described by the following equation. 

Ref .Model - a (standard FEM) + (1-a) (corrected model) 

Thus we can vary the amount of error between the reference model and 
the actual model and monitor the control effort and system response as 

the error increases (a = 0 Perfect Modeling , a = 1 Max. Error) . 

Damping ratios for the simulated plant were obtained experimentally, 
while the damping ratios for the reference model specify the desired 
performance (see table 3) . 

In addition to varying the parameter a, we also varied the initial 
conditions of the structure. In the first simulation the structure 
was deformed into the second mode shape (see figure 6) and released. 
Table 4 shows that introducing errors into the reference model had a 
significant effect on the maximum force at joint 3 and the quadratic 
control effort. The maximum force required at joint 3 is 10 times 
larger at a = 1.0 than at a = 0.0, and the quadratic control effort is 
increased by a factor of 27 over the same interval. In the second 
simulation the structure was deformed into the theoretical fourth mode 
shape (see figure 8) and released. Table 5 shows that introducing 
errors in this case also causes dramatic increases in control costs. 
For example, at the point of maximum error (a = 1), the total control 
effort needed increased by almost 2500%, the maximum control force 
required by actuator 3 increased over 600%, and the amount of^ 
potential and kinetic energy in the system increased over 550%. 


CONCLUDING REMARKS 

The performance of Model Reference Adaptive Control (MRAC) was 
studied in numerical simulations with the objective of understanding 
the effects of differences between the plant and the reference model. 
MRAC was applied to two structural systems with controlled error 
between the reference model and the actual plant. Performance indices 
relating to control effort and response characteristics were monitored 
in order to determine what effects small errors have on the control 
effort and performance of the two systems. It was shown that 
reasonable amounts of error in the reference model can cause dramatic 
increases in both the control effort and and response magnitude (as 
measured by energy integrals) of the plant. 
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TABLE 1 


SUMMARY OF PERFORMANCE MEASURES FOR VARYING 
AMOUNTS OF MASS ERROR IN THE SIMPLY SUPPORTED BEAM EXAMPLE 


% ERROR 

MAXIMUM FORCE 

TOTAL QUADRATIC 
CONTROL EFFORT 

0.0 

380 

11 , 422 

10.0 

540 

27,486 

20. 0 

920 

12 , 263 


TABLE 2 

SUMMARY OF THEORETICAL VS EXPERIMENTAL FREQUENCIES 
FOR THE SLEWING GRID STRUCTURE 


MODE # 

FREQUENCIES (Hz) 


THEORETICAL 

EXPERIMENTAL 

% ERROR 

1 

0.36 

0.42 

16.67 

2 

1.37 

1.45 

5.84 

3 

3 . 00 

2.88 

4 . 00 

4 

4.47 

5.39 

20 . 58 

5 

6.02 

6.41 

6.48 

6 

6 . 69 

6.88 

2.84 | 

7 

9.79 

9.05 

7 . 56 

8 

11.52 

10.18 

11.63 

9 

13 . 11 

13.56 

3.43 

10 

15.35 

14.90 

2.93 

11 

21.16 

15.37 

27.36 
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TABLE 3 


DAMPING RATIOS FOR THE MODELS OF THE SLEWING GRID 


MODE # 

EXPERIMENTAL 

FREQUENCY 

(HZ) 

EXPERIMENTAL 
DAMPING RATIOS 

DAMPING RATIOS 
FOR THE 
REF. MODEL 

1 

0.42 

0.110 

0. 15 

2 

1.45 

0.015 

0.05 

3 

2.88 

0.011 

0 . 05 

4 

5.39 

0.008 

0.05 

5 

6.41 

0.003 

0.05 

6 

6.88 

0.011 

0.05 

7 

9.05 

0.003 

0.05 

8 

10. 18 

0.003 

0.05 

9 

13 . 56 

0.002 

0.01 

10 

14 . 90 

0.002 

0.01 

11 

15.37 

0.002 

0. 01 


TABLE 4 


SUMMARY OF PERFORMANCE MEASURES FOR VARYING AMOUNTS 
OF ERROR IN THE SLEWING GRID 
INITIAL CONDITIONS = MODE SHAPE 2 


ALPHA 

MAXIMUM FORCE (LBS) 

QUADRATIC 

CONTROL 

EFFORT 

lbs 2 -sec 

KINETIC 

ENERGY 

INTEGRAL 

LB-IN-SEC 

POTENTIAL 

ENERGY 

INTEGRAL 

LB-IN-SEC 

ima 

HI 

JT. 5 

0.00 

0.0182 

0.022 

0.021 

0.00044 

0.116 

0 . 117 

0.25 

0.0596 

0.028 

0.025 

0.00107 

0.114 

0 . 119 

o 

• 

(ji 

o 

0.1090 

0.038 

0.035 

0.00299 

0. 112 

0. 120 

0.75 

0.1520 

0.048 

0.045 

0.00742 

0. Ill 

0. 122 

o 

o 

r— 1 

0.1825 

0.057 

0.058 

0.01250 

0. 110 

0 . 125 
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TABLE 5 


SUMMARY OF PERFORMANCE MEASURES FOR VARYING AMOUNTS 
OF ERROR IN THE SLEWING GRID 
INITIAL CONDITIONS = MODE SHAPE 4 


ALPHA 

MAXIMUM FORCE (LBS) 

QUADRATIC 

CONTROL 

EFFORT 

LBS 2 -SEC 

KINETIC 

ENERGY 

INTEGRAL 

LB-IN-SEC 

POTENTIAL 

ENERGY 

INTEGRAL 

LB-IN-SEC 

Bfl 

eh 

JT. 5 

0.00 

0.44 

0.68 

0.58 

0.35 

0.60 

0.60 

0.25 

0.58 

0.78 

0.64 

0.52 

0.70 

0.72 

0.50 

1.19 

0.93 

0.72 

1.47 

1.09 

1.13 

0.75 

2.42 

1.04 

0.89 

5.20 

2.49 

2 . 55 

1.00 

3.12 

o 

• 

H 

0.90 

9.05 

3.94 

3.93 


























Reference Model 



Figure 1 

MODEL REFERENCE ADAPTIVE CONTROL BLOCK DIAGRAM 
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VELOCITV SENSOR 
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FORCE ACTUATOR 


Figure 2 

SIMPLY SUPPORTED BEAM WITH 
VARIABLE CONCENTRATED MASS 
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ce Model Response 



VELOCITY RESPONSE AT THE SENSOR LOCATION FOR THE 
BEAM STRUCTURE WITH A 10% CONCENTRATED MASS 





max no 
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Figure 4 

VELOCITY RESPONSE AT THE SENSOR LOCATION FOR THE 
BEAM STRUCTURE WITH A 20% CONCENTRATED MASS 




* Aluminum Beam Members 
Alloy 6061-T6 
Nominal cross-section: 
2" x 1/8" 

(5.08 cm x 0.32 cm) 


** Ball Bearings 
Make : SKF 


Bearing No. 478207-106 
Pillow Block Flange Unit 


No. FYP-106 


(Bearing seals and all 
grease were removed to 
reduce friction.) 


Figure 5 

SLEWING GRID STRUCTURE 
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